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The effects of density dependence of tlie symmetry energy on tfie collective modes and dynamical 
instabilities of cold and warm nuclear and stellar matter are studied in the framework of relativistic 
mean-field hadron models. The existence of the collective isovector and possibly an isoscalar col- 
lective mode above saturation density is discussed. It is shown that soft equations of state do not 
allow for a high density isoscalar collective mode, however, if the symmetry energy is hard enough 
an isovector mode will not disappear at high densities. The crust-core transition density and pres- 
sure are obtained as a function of temperature for /^-equilibrium matter with and without neutrino 
trapping. An estimation of the size of the clusters formed in the non-homogeneous phase as well as 
the corresponding growth rates and distillation effect is made. It is shown that cluster sizes increase 
with temperature, that the distillation effect close to the inner edge of the crust-core transition is 
very sensitive to the symmetry energy, and that, within a dynamical instability calculation, the 
pasta phase exists in warm compact stars up to 10 - 12 MeV. 
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I. INTRODUCTION 

Understanding the properties of isospin rich nuclear 
matter is one of the present goals of both nuclear physics 
and astrophysics. A major scientific effort is being carried 
out at an international level to study experimentally the 
properties of asymmetric nuclear systems and probe the 
behavior of the symmetry energy close to and above satu- 
ration density Jj . Astrophysical observations of compact 
objects are also a window into both the bulk and the mi- 
croscopic properties of nuclear matter at extreme isospin 
asymmetries [^. Measurements coming from both labo- 
ratory and astrophysical observations are expected to put 
constraints on the acceptable properties of the equation 
of state (EOS) of asymmetric nuclear matter [J, |j| • 

In the last years an important effort has been done 
to determine the density dependence of the symmetry 
energy of asymmetric nuclear matter (see the reviews 
Ul, 0, 0] and references there in) . Correlations between 
different quantities in bulk matter and finite nuclei have 
been established. Examples are the correlation between 
the slope of the pressure of neutron matter at p = 0.1 
fm"'^ and the neutron skin thickness of ^"^Pb ^5], or 
the correlation between the crust-core transition den- 
sity and the neutron skin thickness of ^°*Pb lq|. Both 
Skyrme forces and relativistic mean-field (RMF) mod- 
els have been used to describe these correlations. Re- 
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of the different models and a thermodynamic approach 
to determine the transition density. Many other works 
have studied the instability of asymmetric matter using 
a thermodynamical approach [10|. 

Instead of using a thermodynamic approach, in jll| 
the transition density was estimated in a local equilib- 
rium approximation by determining the density for which 
matter becomes unstable to small density fiuctuations. 
It has been shown in |12h15| that the thermodynamical 
method gives a good estimation of the transition den- 
sity for cold /3-equilibrium matter, although a bit too 
large. In [ij, |lj| a comparison was done between the 
transition density obtained from the Vlasov formalism, 
a semi-classical limit of the description of the collision- 
less regime, and a Thomas-Fermi calculation of the pasta 
phase. It was shown that both results almost coincide: 
the Vlasov formalism predicts a transition density only 
slightly smaller then the Thomas-Fermi one. Moreover, 
Ducoin et al. [lg| have shown that the transition den- 
sity obtained within the Vlasov formalism does not differ 
from the one obtained by [llj using a local equilibrium 
approximation. 

The main goal of this work is to study, in the frame- 
work of relativistic models, the effect of the symmetry en- 
ergy density dependence on some dynamical properties of 
asymmetric nuclear matter, namely the dynamical insta- 
bilities at sub-saturation densities at zero and finite tem- 
oerature and the isoscalar and isovector collective modes 



As shown in [1^ [20| the pressure at the inner bound- 
ary of the crust defines the mass and moment of inertia 
of the crust. This estabhshes a relation between the EOS 
and compact stars observables. Until now the transition 
pressure and density have only been studied at zero tem- 
perature but it is also of interest to determine how do 
these quantities depend on the temperature. This will 
be one of the aims of the present work. The transition 
density and pressure will be determined from the cross- 
ing of the /3-equilibrium EOS of stellar matter and the 
dynamical spinodal. Warm stellar matter may contain 
trapped neutrinos [21|. We will consider the conditions 
expected in the inner crust and in the core of a neutron 
star, namely of isospin asymmetry, density and temper- 
ature. This same problem was investigated very recently 
within a momentum-dependent effective interaction for 
nucleons |22|. 

Collective modes in asymmetric nuclear matter have 
already been studied in previous works J23|-[2»| both 
within non-relativistic and relativistic nuclear models. 
These modes corresponding to isospin or density waves 
in infinitely extended nuclear matter can be identified 
as the counterpart of giant resonances in the atomic nu- 
cleus ^31 ■ They may be of particular importance for the 
formation of shock waves in nuclear matter [2J, [23| and 
they will affect the thermal properties of neutron stars, 
since the neutrino mean free paths are affected by the 
propagation of collective modes in the medium [30j . 

The authors of Ref.[26] have predicted the existence 
of a neutron wave at reasonably large densities, for p ~ 
3poi where po is the saturation density. They have also 
shown that the inclusion of the scalar isovector virtual 
(5(ao(980)) field may strongly reduce the restoring force 
for the isovector collective modes of the system. In a 
similar way they predict that, if the incompressibility is 
decreasing with the increase of the density, the exotic 
high density isoscalar collective modes could disappear. 
This could be a signature of the softening of the equation 
of state (EOS) at large densities [2g . 

In the present work we will study the collective modes 
of nuclear matter within the Vlasov equation for rela- 
tivistic nuclear models delevoped in [3l|, [33] for symmet- 
ric nuclear matter and applied later to asymmetric nu- 
clear matter [23, [2^ |3j] . It is a semi-classical approach 
which can be used to determine the eigenmodes of nu- 
clear asymmetric matter and corresponds to the fi ^^ Q 
limit of the time dependent Hartree-Fock. The relativis- 
tic Vlasov equation for the non-linear Walecka model 



that differ in the density dependence of the symmetry 
energy and the incompressibility. The following models 
will be considered: NL3 [3a] with a quite large sym- 
metry energy and incompressibility at saturation and 
which was fitted in order to reproduce the ground state 
properties of both stable and unstable nuclei, TMl [S^ 
which also reproduces the ground state properties of both 
stable and unstable nuclei and provides an equation of 
state of nuclear matter similar to the one obtained in 
the RBHF (Relativistic Brueckner Hartree-Fock) theory, 
softer than NL3 at high densities, FSU [3^ which was ac- 
curately calibrated to simultaneously describe the GMR 
in 90Zr and ^ospb, and the IVGDR in ^ospb and still 
reproduce ground-state observables of stable and unsta- 
ble nuclei. In TMl a quartic self-interaction w-meson 
term was included, which softens the EOS. This term is 
also present in the FSU parametrization together with 
a mixed isoscalar-isovector coupling which modifies the 
density dependence of the symmetry energy. Finally, we 
also consider NL3wp and TMlwp, the NL3 and TMl 
parametrizations with a mixed isoscalar-isovector cou- 
pling which we will vary in order to change the density 
dependence of the symmetry energy [^. The EOS ob- 
tained within NL3 is too hard and well above the lim- 
its that collective flow data [S^] impose to the EOS of 
symmetric matter. Also, the symmetry energy slope of 
NL3 at saturation is 118 MeV slightly above the limit 
obtained with experimental constraints from isospin dif- 
fusion 63 < L < 113 MeV, and well above other con- 
straints [i^. Nevertheless, we have chosen this model 
because it has been calibrated to describe the ground- 
state of a wide number of stable and unstable nuclei and 
because it is still frequently used. This model is taken 
as an example of a model with an hard EOS both in the 
isoscalar and isovector channel. 

We will consider both asymmetric nuclear matter, neu- 
tral neutron-proton-electron [npe) neutrino-free matter 
in /3-equilibrium at zero temperature and npe matter 
with trapped neutrinos in /3-equilibrium for a lepton frac- 
tion Yi = 0.4, both at zero and finite temperature. 

In Sec. |TT] we review the formalism used, in Sec. IIIII 
we present and discuss the results obtained and finally, 
in Sec. IIVI we draw some conclusions. 



II. THE VLASOV EQUATION FORMALISM 



We use the relativistic non-linear Walecka model 



the electromagnetic field A^. The Lagrangian density 
reads: 

where the nucleon Lagrangian reads 

with 
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and the electron Lagrangian is given by 

Ce = i'e [it, {id^ + e^^) - me] Ve. 

The isoscalar part is associated with the scalar sigma 
{(j) field (/), and the vector omega (w) field V^, whereas 
the isospin dependence comes from the isovector-vector 
rho {p) field 6* (where /i stands for the four dimensional 
space-time indices and i the three-dimensional isospin 
direction index). The associated Lagrangians are: 
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where O^j, = 5^K - d^V^,, B^^, = d^h^ - d^h^, - 
gpihfj, X h^) and F^,^ = df,A^ - d^A^. 

The model comprises the following parameters: three 
coupling constants gs , gv and gp of the mesons to the nu- 
cleons, the bare nucleon mass M, the electron mass me, 
the masses of the mesons, the electromagnetic coupling 
constant e — ■y/47r/137 and the self-interacting coupling 
constants k, A and ^. In this Lagrangian density, t is the 
isospin operator. 

The NLSwp, TMlwp and FSU parametrizations also 
include a nonlinear wp coupling term, which allows to 
change the density dependence of the symmetry energy. 
This term is given by: 



Table |T] shows nuclear matter properties at the sat- 
uration density for these models: the binding energy 
per nucleon E/A, the incompressibility coefficient K, 
the symmetry energy Ssym, the symmetry energy slope 
L, the symmetry energy curvature Ksym and AV ~ 
Ksym -6L-Q^L (see Ref.jll]). 



We have chosen a set of models which have different 
incompressibilities at large densities and different density 
dependences of the symmetry energy. This is clearly seen 
in Fig. [T] where we have plotted the pressure of symmetric 
nuclear matter, Fig. [T^), and the symmetry energy for 
densities below 0.45 fm~^. Fig. [TJj). In Fig. [Iji) we also 
include the constraints obtained from the flow analysis in 
heavy ion collisions [39]. It is seen that NL3 and NL3wp 
models give a too hard EOS at large densities, which 
is not compatible with the constraints from heavy ion 
collisions. FSU is well inside the constraints and TMl 
and TMlwp are just slightly above the upper limit. 



The non-linear cup term will only affect the symme- 
try energy, therefore, models that only differ by an extra 
non-linear top term will have the same EOS for symmetric 
matter. As we will show later, this behavior will have a 
strong effect on the high energy collective modes. With 
respect to the symmetry energy NL3 and TMl have a 
very hard behavior at high densities while FSU is softer. 
The NL3ujp and TMluip models show an intermediate 
behavior and approach FSU when the cop coupling pa- 
rameter A^ increases. Below p = 0.1 fm~^ the softness 
of the symmetry energy is reversed with NL3 the softest 
and FSU the hardest. 



In the sequel we use the formalism developped in ref. 
[32l | where the collective modes in hot and dense sym- 
metric nuclear matter were determined within the Vlasov 
equation based on the Walecka model [44|. We will use, 
whenever possible, the same notation. We start from the 
thermodynamical potential density for npe matter 
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FSU 




0.148 


-16.302 


227.895 32.537 


60.395 


-51.415 


-275.570 


TMl 




0.145 


-16.265 


279.514 36.835 


110.607 


33.557 


-517.390 


TMla;p,A„ = 


= 0.01 


0.145 


-16.265 


279.514 34.855 


85.290 


-79.145 


-503.986 


TMla;p,A„ = 


= 0.02 


0.145 


-16.265 


279.514 33.235 


67.654 


-97.150 


-434.142 


TMla;p,A„ = 


= 0.03 


0.145 


-16.265 


279.514 31.880 


55.082 


-72.164 


-346.532 
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0.148 


-16.240 


269.937 37.344 


118.320 


100.525 


-696.129 


NL3ajp,A„ = 


0.01 


0.148 


-16.240 


269.937 34.929 


87.638 


-46.2956 


-636.362 


NL3ajp,At, = 


0.02 


0.148 


-16.240 


269.937 33.121 


68.147 


-53.458 


-512.296 


NL3ajp,At, = 


0.03 


0.148 


-16.240 


269.937 31.660 


55.227 


-8.069 


-379.915 



TABLE I. Nuclear matter properties at saturation density, po- AH quantities are in MeV, except for po, given in fm 
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FIG. 1. (Color online) a) Pressure of symmetric nuclear matter for NL3, TMl and FSU and b) symmetry energy for NL3 and 
NL3a;p, TMl and TMlup, and FSU. We have considered A„ = 0.01, 0.03 in NL3a;p and TMlojp. For symmetric matter the 
EOS of NL3ajp and TMlup coincide respectively with NL3 and TMl. In a) we also include the constraints obtained from the 
flow analysis in heavy ion collisions (shaded region) j3^] . 



where £ is the energy density 
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to (Vi, 6i, Ai ), <; is the entropy density 
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and /ii and pi are the chemical potential and particle 
density of particle type i. We only consider the third a^ 
component in isospin space of the p- meson and we denote 

by 

ffr,n,i)+=diag(f.+ . f„+, f.+ ) 



For protons and neutrons, i = p,n, we have 



where 



Voi = gvVo + Y'^ibo + eAo — — ^, 

with Ti — 1 (protons) or —1 (neutrons). For electrons, 
i — e, we have 
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The time evolution of the distribution functions is de- 
scribed by the Vlasov equation 



dt 



+ {/»±>»±} = 0, 



:p, n, e. 



(7) 



where {, } denote the Poisson brackets. The Vlasov equa- 
tion expresses the conservation of the number of particles 
in phase space and, therefore, at a given temperature, 
the dynamics described by the Vlasov equation does not 
change the number of particles or antiparticles in this 
semiclassical approach. 

From Hamilton's equations we derive the equations de- 
scribing the time evolution of the fields (/>, V^ , A^^ , and 
the third component of the p- field, 63 = (&01 b), which are 
given in Appendix El (Eqs. (|M|) - (|X71) V The state which 
minimizes the thermodynamical potential density of npe 
matter is characterized by the distribution functions 
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The chemical equilibrium conditions /ip = /i„ — /-ie for 
neutrino free matter, or /Xp = /x„ — ^g -f ^.^^ for neutrino 
trapped matter, as well as the electrical charge neutrality 
condition pe = Pp must be imposed. 

Collective modes in the present approach correspond 
to small oscillations around the equilibrium state. These 
small deviations are described by the linearized equations 
of motion and collective modes are given as solutions of 
those equations. To construct them, let us define: 

f^±=f^■l+Sf^±, 

bo = b^o"' +Sbo, bi = Sbi, 
Ao = SAo, A, = SA,. 
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and by the constant mesonic fields, which obey the fol- 
lowing equations: to^^q + f'/'o + ^Vo = dsPs , 
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As in [27], |31|, [32, |4J| j we express the fluctuations of the 
distribution functions in terms of the generating func- 
tions: 



S±{r,p,t) = diag(S'p±, S„±, Se±) 



such that 



6f± = {S±Jo,±}^{S±,p^} 



dfo± 
dp^ 



The linearized Vlasov equations for 5fi± , 
d5fi± 



dt 



-l-{(5/i±,/ioi±} + {/oi±,(5/ii±} = 



are equivalent to the following time-evolution equations: 



where 

^Voi = QvSVq + Ti ^ (56o + e — — ^ Mq, 

with hm± = ±eoi + Vg^' and /ioe± = ±eoe- 

Of particular interest on account of their physical rel- 
evance are the longitudinal modes, with wave vector k 
and frequency w, described by the ansatz 
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The eigenmodes uj of the system, with sound velocity 
Vs = w/fc, are solutions of the Eq. (IT71) . At subsatu- 
ration densities, there are unstable modes identified by 
imaginary frequencies. For these modes we define the 
growth rate T — —icu. The region in {pp,Pn) space for a 
given wave vector k and temperature T, limited by the 
surface uj — 0, defines the dynamical spinodal surface. 
In the fc = MeV limit we recover the thermodynamic 
spinodal which is defined by the surface in the {pp, p„, T) 
space for which the curvature matrix of the free energy 
density is zero, i.e., has a zero eigenvalue. This relation 
has been discussed in 1331. 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 



where j = p, n, e, C = V, i>, A represent the vector-meson 
fields and 6 is the angle between p and k. The wave 
vector of the excitation mode, k, is identified with the 
momentum transferred to the system through the process 
which gives rise to the excitation. 

For the longitudinal modes, we get 5V^ — SV^ — , 
5bZ = Sbl = and SAZ ^ SAv = . Calling 6V/: = 
SVuj, Sb^ — Sb^ and SA^ — SA^, we will have for the 

fields 6V^,^ = JV«_e'("t-kr) ^^^^ gy^^ ^ JyMg»(a;t-k.r)^ 

Replacing the ansatz (TTSl) in Eqs. (fTTj) and (|T2|) we get 
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In Appendix |B] we present the equations of motion for 
the field fluctuations (Eqs. ((BT|) - ((B7)) '). The solutions of 
these equations and Eqs. (|14p . ([15]) form a complete set 
of eigenmodes that may be used to construct a general 
solution for an arbitrary longitudinal perturbation. Sub- 
stituting the set of equations (|BT]) - ((B41) into dUl) and p^ 
we get a set of equations for the unknowns 5^^, which 
lead to the following matrix equation 
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In the present section we discuss the effect of the sym- 
metry energy and incompressibility on some of the dy- 
namical properties of nuclear and stellar matter. We cal- 
culate the collective modes at densities above and below 
saturation density and analyze their isoscalar or isovector 
character. We will next determine the dynamical spin- 
odal surfaces, for different temperatures and wave vec- 
tors, characterized by a zero frequency solution of the 
dispersion relation (1171) . They will allow us to predict 
under which conditions a non-homogeneous phase at the 
crust of a proto-neutron star may exist. We will finally 
determine for warm /3-equilibrium matter with trapped 
neutrinos the unstable modes with the largest growth 
rate, which define the mode that drives the system to 
a non-homogeneous phase and, therefore, gives an esti- 
mation of the size of the clusters formed in the phase 
transition. For these modes the magnitude of the distil- 
lation effect will be calculated. For /3- equilibrium warm 
stellar matter with trapped neutrinos we will consider a 
fixed lepton fraction Yi = 0.4, which corresponds to the 
maximum expected value [2ll |45| . 



A. Collective modes 

We have calculated at zero temperature the normal 
modes of asymmetric nuclear matter as a function of the 
nuclear matter density p and for different proton fractions 
Up = Pp/P- Fig. [2] shows in the top row the ratio of the 
proton to neutron density fluctuations, Eq. (jC3|) . and 
in the bottom row the collective mode sound velocitv. 
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FIG. 2. (Color online) The ratio of the proton to neutron density fluctuations [Eq. (|C3|) ] divided by the ratio of the proton 
to neutron densities [(a) j/p = 0.1, (b) yp = 0.3 and (c) j/p — 0.4] and the sound velocity Vs normalized to the neutron Fermi 
velocity, s„ = Vs/vFn, [(d) yp = 0.1, (e) yp = 0.3 and (f) yp = 0.4] for NL3 and NL3a;p, TMl and TMlcjp, and FSU. A„ = 0.01 
and 0.03 have been used in NL3ajp and TMlujp. The curves have been obtained for T = MeV and k = 200 MeV. 
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FIG. 3. (Color online) Isovector-like mode below 2po for 



above twice the saturation density, which we will des- 
ignate by crossing density Pc- The crossing density de- 
creases with an increase of the proton fraction. This is 
clearly seen in Fig. |4]where the crossing density is plotted 
for NL3 and two parametrizations of NL3wp and TMlwp 
(A« — 0.01 and 0.03). The crossing density is sensitive 
to the high density dependence of the symmetry energy 
and, therefore, the larger the crossing density the larger 
the difference between NL3 and NLSwp. We point out 
that NL3 and NL3a;yO only differ on the isovector channel 
and therefore, the differences among these models reflect 
the isovector behavior. TMl has, like NL3, a very stiff 
symmetry energy, but, unlike NL3, has a quite soft EOS 
for symmetric matter at large densities due to the inclu- 
sion of quartic w-terms. As a result the collective mode 
exists at large densities but with an isovector character. 
Including the non-linear top terms in TMl changes this 
situation: for a strong enough coupling the isovector col- 
lective mode will change into an isoscalar mode, however, 
at much larger densities than NL3 and NL3a;p models, 
see Fig. |4l 

At very low densities the sound velocity of the modes 
decreases as the asymmetry increases and, as seen in 
FiglSk) for 2/p = 0.1 there are no undamped modes for 
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FIG. 4. (Color online) The crossing density pc defining the 
onset of isoscalar character in the collective modes in np mat- 
ter versus the proton fraction, for k=200 MeV and the models 
NL3, NL3a;p and TMla;p with A„ = 0.01, 0.03. 
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high densities due to the inclusion of the (S-meson. In a 
similar way the isoscalar mode depends on the isoscalar 
response of the system and may not exist if the EOS is 
too soft. 

The properties of the collective modes of nuclear mat- 
ter described by the models introduced in the first two 
sections are determined by the interaction factors fi 

\ 2 

9i 



h 



s, V, p, 



where the effective meson masses mi^eff are given by 



(18) 



's.e// 



ml + k4)o + -4>l + g 



f dps 



V dM* 



2 2,4 4ta2 I f) a 2 2i2 

2„2ta2 



(19) 

(20) 

(21) 

At T = MeV, the dispersion relation pT)) may be 
reduced to a simple analytical form for symmetric matter 
and in the low momentum limit |26| . If we make the 
approximation Vg ~ vp — kp/^F we get for the isoscalar 
mode 



eff = K + 2A.r7'5p^o • 



1 + 3|^P 



M* 



ifv - fs) (1 - Vp) + fyfs Vp^P, 



and for the isovector mode 






0s =0 

(22) 
(23) 



where (ps is the Lindhard function (ps ^ 1 — s/2 In |(s + 

1)/(S- l)l> S = Vs/Vp. 

We notice that the second term of Eq. ([2^ is dom- 
inated by the first part propotional to {f^ — fs), while 
the second term of Eq. (|^5|) strongly depends on fp. 
Both these factors are multiplied hy 1 — Vp — {M* /epY 
which approaches fast zero when the density increases. 
In Fig. [5] we represent these quantities: (/.„ — fs) and 
{fv - fs){l - v\) in Fig. EJi) and fp and fp{l - vj,) in 
Fig. Eb) for NL3, TMl and FSU, and /p(l - vj,) for 
TMlcjp and NLSwp in Fig. [St). There will be a isoscalar 
mode only if the factor that multiplies the Lindhard func- 
tion is positive. This only happens when the repulsion 
described by the oj-meson overcomes the cr-meson attrac- 
tion, which occurs for p > 0.2 fm~'^ for NL3 and TMl 
and for p > 0.3 fm~'^ for FSU. However, the increase of 
the effective mass of the w-meson in TMl and FSU due 
to the quartic w-term in the Lagrangian density [see Eq. 



and the response of the system to isovector excitations 
([23|) . which depend on fp, wiU be significantly softened 
at large densities. Moreover, since the gp parameter is 
fixed so that at p = O.f fm^'^ the symmetry energy is 
the same for the models with and without the non-linear 
cop term, the larger the coupling parameter A^ the larger 
must he gp. A larger gp enhances the value of fp at small 
densities, and, therefore, also the response of the system 
to isovector perturbations as seen in Fig. |3j the models 
with the larger sound velocity and the density fluctuation 
ratio close to -1 are the ones with the larger fp. 

Comparing fp with fp{l—Vp) [Fig. [5}))], it is clear that 
the factor l~Vp makes the second term of pB]) go faster 
to zero at large densities. TMl has the largest values 
of /p(l — w|,) at large densities which may explain its 
different behavior at these densities, namely the presence 
of an isovector like mode. 

We conclude that the main role of including the non- 
linear tup term in NL3 and TMl is to reduce the value 
of fp at large densities, making the isovector character of 
the modes weaker in this density range, and to increase 
the fp at low densities giving rise to stronger isovector 
modes below p « 0.2 fm~^. 

We stress that the above discussion refers to symmet- 
ric matter at zero temperature. In asymmetric matter 
isoscalar and isovector modes are coupled and we ex- 
pect to find the behavior of both pure isoscalar and pure 
isovector modes in all modes. A stronger fp (/^ — fs) 
will give a stronger isovector (isoscalar) character to the 
mode. 

In Fig. [HI we present the proton fractions (a) the 
proton-neutron density fluctuations (b) and the corre- 
sponding sound velocities of the collective modes s„ — 
Vs/vpn (c) for /3-equilibrium cold stellar matter. For 
the FSU parametrization we get no modes and for the 
parametrizations TMl and TMlojp only isovector-like 
modes are present. The other models change charac- 
ter from isovector to isoscalar at the crossing density pc. 
The crossing density in /3-equilibrium matter is almost 
the same for all NL3 like models although, it was seen in 
Fig. [2] that the crossing density depends on the isovec- 
tor properties of the model. This is due to the fact that 
the models have different proton fractions for the same 
density and there are two opposite effects that define 
the crossing density: the larger the proton fraction the 
smaller the crossing density and the larger the symmetry 
energy at saturation the larger the crossing density. So 
the combination of these two effects leads to a similar 
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FIG. 6. (Color online) /3-equilibrium stellar matter: (a) pro- 
ton fraction (b) the ratio of the proton-neutron density fluc- 
tuations (Eq. (|C3|l ) and (c) the corresponding neutron sound 
velocities s„ for NL3, FSU, TMl, and NL3a;p and TMlup 
with A^ = 0.01 and 0.03. The results are for T = MeV and 
k = 200 MeV. 



odal. The extension of the spinodal region depends on 
the wavelength of the perturbations: the system is sta- 
ble to small wavelength perturbation, of the order of the 
range of the nuclear force or smaller and to large wave- 
lengths due to the contribution of the Coulomb interac- 
tion which blows up at large wavelengths. Fig. [7] shows 
the approximate envelope of the spinodal regions for dif- 
ferent temperatures and for NL3, FSU, TMl and NLScjp 
with A^ = 0.01,0.02,0.03. We will not consider TMl 
with non-linear ujp because at subsaturation densities 
TMl and NL3 behave in similar ways. With the increase 
of temnerature. the sninodal reffion shrinks to a, reffion 



the symmetry energy. 

The transition densities at the crust-core transition for 
/3-equihbrium warm stellar matter with trapped neutri- 
nos (Yi = 0.4) are shown in Table HIl for several temper- 
atures. These densities are obtained from the crossing 
of the EOS of /3-equilibriuni matter with the dynamical 
spinodal. The transition densities decrease with tem- 
perature due to the decrease of the spinodal region. At 
T —0 MeV, there is no clear dependence of the transi- 
tion densities on the slope L. However, if we consider 
the set of models that has the same isoscalar description, 
NL3 and NLSojp, a clear tendency is seen: pt decreases 
when L increases. This is due to the corresponding de- 
crease of the transition proton fraction (for T —0 MeV 
yp = 0.312, 0.314, 0.318, 0.322, respectively, for NL3, 
NL3wp with A^ = 0.01, 0.02,1,0.03). A larger L cor- 
responds to a smaller symmetry energy at subsaturation 
densities and allows a larger asymmetry. 

In Table nil we also show the transition densities for [3- 
equilibrium neutrino-free stellar matter at T = MeV. 
For neutrino-free matter, the transition densities are 
smaller than the ones with trapped neutrinos because 
neutrino free matter has a much smaller proton fraction. 

The pressure has no monotonic behavior with temper- 
ature: below T ^ 5 MeV it increases with temperature 
and only above T ^ 5 — 10 MeV it reduces with temper- 
ature due to a large reduction of the transition density. 
At lower temperatures the small reduction of the tran- 
sition density with respect to the T = MeV case and 
a larger lepton contribution to the pressure explains the 
increase of the transition pressure. Similar results have 
been obtained in }22] . 

The effect of the symmetry energy slope L on the tran- 
sition pressure of matter with trapped neutrinos is only 
clear for the set of models with the same isoscalar EOS, 
namely NL3 and NLSojp. For these models the transition 
pressure increases with L at T =0 MeV but does not show 
a systematic behavior at T =5 and decreases with L at 
T = 10 MeV. At T =0 MeV, the smaller isospin asymme- 
try for the smaller L corresponds to a smaller pressure. 
This tendency inverts as the temperature increases due 
to the larger lepton contribution to the pressure for the 
larger densities. 

For neutrino-free matter at T =0 MeV it is seen that 
taking the whole set of models the transition pressure 
does not show a clear dependence with density, in agree- 
ment with the findings of [46] with a much larger set of 
models. However, if we limit ourselves to the set NL3 and 



pointed out that the present method of estimating the 
cluster size takes into account both the Coulomb inter- 
action, which has been included in the formalism, and 
surface effects, which are due to the finite range of the 
nuclear force. Within a relativistic mean-field approach 
to nuclei and finite clusters surface properties are deter- 
mined by the exchange of mesons with a finite mass. The 
estimation done with the present method agrees with the 
results obtained within a Thomas Fermi calculation of 
the pasta phase [l3|. 

In Fig. |5]we show the cluster size estimates for matter 
in /3-equilibrium at zero temperature. Large differences 
for different models can be seen for the larger densities. 
This is due to the fact that cold /?-equilibrium matter 
has quite small proton fractions, for which larger differ- 
ences arise between the models. It is also interesting to 
notice that close to the crust-core transition the clus- 
ters get larger and all models predict a size of the order 
of 10 fm. There is also a close correlation between the 
extension of the non-homogeneous phase determined by 
the density at which each curve ceases to exist, and the 
slope of the symmetry energy as noticed by many authors 
7, 48]. While below p = 0.02 fm"^ all the models pre- 
dict the same size of clusters, above this density there are 
differences between the models which are reflected in the 
fact that all predict similar sizes at the transition density, 
although this density differs from model to model. 

In Fig. iniwe plot the proton over neutron density fluc- 
tuations, normalized to their corresponding Pp/pn ratios 
for the largest unstable mode (flrst row) , the correspond- 
ing growth rates (second row) and the estimated size 
of the clusters (third row) in /3- equilibrium stellar mat- 
ter with trapped neutrinos for T = 5, 10 and 12 MeV. 
These unstable isoscalar modes, since protons and neu- 
trons fluctuations are in phase, are the ones that drive the 
system to a separation into a dense symmetric phase and 
a neutron rich gas. They are different from the isovec- 
tor modes represented in Fig. [2] at subsaturation densi- 
ties, for which neutron fluctuations are larger than proton 
fluctuations. This property has already been point out 
in [2^. 

When the temperature is increased, the ratio of proton- 
neutron density fluctuations decreases as well as the den- 
sity range for clusterization. On the other hand, there 
is an increase of the average size of the clusters. From 
the top row of Fig. [HI we conclude that for warm stellar 
matter with trapped neutrinos, FSU and NLSojp are less 
effective in recovering; the isospin svmmetrv in the liquid 
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FIG. 7. (Color online) Spinodals (Eq.dTTll, with cj = 0) for fe = 75 MeV for NL3, FSU, TMl and NL3a;p with A„ = 0.01, 0.02 
and 0.03 at a) T = 0, b) T = 10 and c) T = 12 MeV. 



r(MeV) 



Pt (fm-^) 



Pt (MeV/fm'^) 



NL3a;p 
A„ 
FSU TMl NL3 0.01 0.02 0.03 FSU TMl NL3 



0.01 



NLSwp 
A„ 
0.02 0.03 



y„=0 0.0745 0.0602 0.0538 0.0619 0.0741 0.0855 0.388 0.324 0.236 0.357 0.500 0.535 



y;=0.4 0.0814 O.O829 O.O8I6 O.O828 0.0840 0.0850 0.665 0.844 0.865 

5 0.0789 0.0803 0.0781 0.0795 0.0809 0.0821 1.021 1.064 1.087 

10 0.0652 0.0671 0.0605 0.0621 0.0642 0.0667 1.059 1.097 0.971 

12 0.0475 0.0545 - - _ _ 0.842 0.969 - 



0.838 0.802 0.732 
1.074 1.086 1.090 

1.022 1.061 1.102 



TABLE II. The transition densities and the pressures for NL3, TMl, FSU and NL3a;p with A„ = 0.01, 0.02, 0.03 at the crust- 
core transition for /^-equilibrium neutrino-free stellar matter (Y^=0) at T = MeV and for /3-equilibrium stellar matter with 
trapped neutrinos {Yi — 0.4) for several temperatures. 
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[l3j). There, a simple approach was used to determine 
the clusters: a zero surface thickness ansatz for the clus- 
ters was used with a non self-consistent surface energy. 
It is, therefore, important to understand this point by 
describing the clusterized phase within a finite temper- 
ature self-consistent Thomas-Fermi calculation. In the 
approach used here the increase of cluster sizes with tem- 
perature is a consequence of the fact that the smaller 
wavelengths do not give rise to instabilities when tem- 
perature increases. At T = 12 MeV only TMl and FSU 
models predict the formation of clusters. Globally, the ef- 
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FIG. 9. (Color online) /3-equilibriuni neutrino-trapped matter {Yi = 0.4): the normalized ratio of the proton over the neutron 
density fluctuations for the most unstable mode (largest |aj|) at (a) T = 5 , (b) T = 10 and (c) T — 12 MeV; the corresponding 
growth rates at (d) T = 5 , (e) T = 10 and (f) T = 12 MeV and estimated sizes of the clusters at (g) T = 5, (h) T = 10 and 
(i) T = 12 MeV for NL3, FSU, TMl and NL3a;p with A^ = 0.01, 0.02 and 0.03. 




P( fm-^) 



d£. 



FIG. 10. (Color online) The ratio £7^—5— 
tion densities for the models under study. 



0.08 



at subsatura- 



where J-^ is the second derivative of the free energy den- 
sity of symmetric matter with respect to the density. 
Since the isoscalar contribution of the different models 
does not differ much at subsaturation densities we may 
expect that the different behavior between the models 
is determined by the ratio p{d£sym,/dp)/£sym- We have 
plotted this quantity for the range of densities of interest 
to understand the nonhomogenous phase. We conclude 
that this quantity seems also to define the behavior of 
the ratio of the proton to neutron density fluctuations of 
the most unstable mode in a dynamical approach to the 
instabilities of asymmetric nuclear matter. 



IV. CONCLUSIONS 



much. 



In the present work we have compared normal modes 
nronerties described bv several rela,tivistic nuclear mod- 



neutral npe matter. We have considered a) collective 
modes at zero temperature, both isovector and isoscalar- 
like, in a wide range of densities and isospin asymmetries 
b) unstable modes at subsaturation densities for both 
cold and warm matter. We have analyzed the effect of 
the isovector properties, namely the density dependence 
of the symmetry energy, on the collective modes. All 
models have an isovector-like mode for densities p < 0.2 
fm~'^ and yp -^ 0.25. However, when the proton frac- 
tion decreases this mode becomes weaker and exists in a 
smaller range of densities. Models with a smaller sym- 
metry energy slope have a stronger mode with a larger 
sound velocity. It was shown that this was due to the 
larger coupling factor fp these models have at subsatu- 
ration densities. 

For p > 0.2 fm~^, FSU has no collective mode. The 
parametrization TMl has a mode which keeps always 
its isovector-like character, also for very high densities. 
For NL3, NLSwp and TMlwp the modes change from 
isovector-like to isoscalar-like for a density above 0.2 
fm^'^ that depends both on the model and the proton 
fraction. At the crossing the collective mode is a neutron 
wave. This feature was first identified in ^26,]. The differ- 
ences among the models are due to the different behaviors 
at large densities of their effective coupling contants fi 
and the nucleon effective mass M* . A finite temperature 
calculation of these modes should be performed because 
it may affect the evolution of a proto-neutron star with 
trapped neutrinos. 

For /3-equilibrium cold neutrino free matter FSU has 
no collective mode, TMl has an isovector mode for all 
densities above ^ 0.11 fm~^, NL3 and NLSwp have a col- 
lective mode respectively above 0.15 and 0.3 fm^^ which 
changes from an isovector to an isoscalar-like mode at 
p ~ 0.4 fm^"^. It is seen that the models considered pre- 
dict quite different dynamical responses. 

Next we have studied the unstable modes which all 
models present at subsaturation densities, both for cold 
and warm matter. The dynamical instability density 
region defines a lower bound for the non- homogeneous 
matter in the crust of a compact star [ij, [Ij, [l^ . We 
have compared the previsions of the different models for 
cold /?-equilibrium matter and warm /3-equilibrium mat- 
ter with trapped neutrinos. It was shown that the non- 
homogeneous phase could extend to ~ 12 MeV in some 
models but at T = 10 MeV all models predicted still a 
pasta phase. While the extension of cold /3-equilibrium 
matter is model denendent and shows a, close correlation 



smaller distillation effect and, therefore, a proton richer 
background gas. This will have an effect at the level of 
transport properties such as electrical conductivity. It 
was also shown that the average cluster size increases 
with temperature. For cold /3-equilibrium matter the av- 
erage cluster size at the crust-core transition density does 
not depend on the model and is 10 fm large. However, be- 
cause the transition density varies with the model there 
is quite a dispersion on the average cluster size above 0.02 
fm~'^: the models with the larger pasta extension predict 
smaller clusters, ~ 6 fm, for a wider range of densities. 
It was also shown that while the transition density de- 
creases systematically with the increase of temperature, 
the transition pressure increases until T ^ 5 MeV and 
decreases for larger temperatures. No systematic behav- 
ior of the transition pressure with the symmetry energy 
slope L was obtained if the whole set of models is consid- 
ered both for neutrino-free matter and neutrino trapped 
matter. However, taking the set of models with the same 
isoscalar EOS, NL3 and NL3wp, we could show that for 
neutrino-trapped matter at T =0 MeV the pressure de- 
creases when L decreases due to the smaller isospin asym- 
metry, at T = 5 MeV there was no systematic behaviour 
due to a competition between the lepton contribution 
and the nucleon EOS and at T = 10 MeV the pressure 
decreases with an increase of L due to the smaller lepton 
contribution, for the smaller densities. 

We conclude that transport coefficients of stellar mat- 
ter will be sensitive to the density dependence of both the 
EOS of symmetric matter and the density dependence of 
the symmetry energy and may be used to set constraints 
on the EOS of cold and warm stellar matter. 



Appendix A: Equations for the fields 
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where the scalar density is 
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the equations of motion for the field fluctuations read 
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The components of the baryonic four-current density are 
d^p 
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and from the continuity equation for the density currents, 
we get for the components of the vector fields 
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Appendix C: Solutions for the eigenmodes and the 
dispersion relation 
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and the components of the isovector four-current density 
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We define the following quantities 
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The coefficients aio are defined as: 



-2 



(27r)2r \k 



IcJe) 



p/e 



dx- 



X 



p/e 



UJ ±X 



± 



2^ 



■ uj In 



u-p/e 



UJ +ple 



jm 



e"/.^(e)/o.T(/o.T-l)&, 



MT 



an = 1 + Cs [ll\ - /°^ j , ai2 
«i3 = -(cr + Ce)(/:^++e) 



7-Op _ tOp 



ai4 = -cr (/ip 



LJ \ o; 



rip 



rip 



-c^- / > ai5 - Ce [I^+ I -.„ 



ip 



a2A^'cl-{lll+llrS), 025-0, 



031 



Ip 



rip 



, ^32 
a33 = l-(c£P + Ce)(/^^+-/, 



rip 



rip 



j2p 



034 



rfe-e),«35 

In ~ 



r2p 



,2p 



a, (Z.^" 



ill 



"43 = -Cj {III ^'-) . 

a44 = 1 - cr {III - C-) , «45 = 0, 
051 — a52 = 054 = 0, 



rln 



s V "+ ' 1^ — / ' 



r2e 



«53 - Ce (e+ - II-) , 055 = 1 - Ce (e+ " I'^-J 

The ratios of the amplitudes are given by: 
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